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\nsr?ACT 


ilcctrochemical machining is one of the most widely 

used unconventional manufacturing process. Its capabilities have 
not been fully exploited mainly due to some inherent problems 
associated with tool design. Conventional IC*' tool design is not 
accurate enough to account for different complex phenomenon taking 
place in interelectroc’o gap. One such phenomenon is. heat transfer 
through the electrodes. There is no accurate model available for 
the hoot transfer analysis of 7Cr. 

In the present work, a two-dimensional, unsteady state 
heat transfer model has been proposed to analyze the heat transfer 
from electrolyte to and through the electrodes, for obtaining the 
temperature distribution in interelec trodo gap. ’Experimental 
data have been used to test the accuracy of the model. Comparison 
of analytical and experimental temoerature distribution reveals a 
good agreement between the two* It has been concluded that the 
anode shape can be predicted more accurately by incorporating 
two-dimensional, unsteady state heat transfer. 

?*.ajority of the models available for tool design in 
ignore the effect of variation in different process parameters. 

A method has been suggested for the determination of the shape of 
the tool (cathode) to produce the required work (anode) shape, 
while machining under specified machining conditions. This method 
takes into account the effect of variation in different process 
parameters. The cathode shape prediction problem has bean attem- 
pted for electrochanical drilling with a bare tool, electrochem- 
ical bit drilling and electrochenical bit boring. Tools and tool 
bits have been designed for producing profiles obtained experi- 
maatally and also for producing straight sided holes* Comparison 
of designed and expeilraental tool shapes reveals a good agreement 
between the two. 



CHAPTER I 


INTRODUCTION AND LITERATURE SURVEY 


1.1. ELECTROCHB^CAL MACHINING (EaOs 

The rapid developments in the field of materials have 
given an impetus to the modem manufacturing technology to 
develop, modify and discover newer technological processes with 
a view to achieve results that are far beyond the scope of 
existing conventional or traditional manufacturing processes. 

As a result of such developments in the field of materials, 
some new techniques of machining such as ECM, EDM, USM, AJM, 
EBM, etc., have been developed. These nontraditional processes 
are capable of providing effective solutions to the problems 
imposed by the increasing demand for accurate machining of 
high-strength-tCTiperature-resistant (HSTR) alloys, the require- 
ments of parts with intricate and complicated shapes, and 
materials so hard as to defy machining by conventional methods. 

Electrochenical machining (ECM) is one of the most 
widely used unconventional machining process. ECM is a process 
of controlled and accelerated anodic dissolution of metals and 
alloys. High velocity electrolyte flows between anode and 
cathode (Figure 1.1) subjected to a small d.c. voltage and 
maintaining a small interelectrode gap (lEG) between them. 

The flowing electrolyte helps not only in allowing the high 
rate of metal dissolution, but also takes away the reaction 
products and heat generated during the chemical and electro- 
ch®fiical reaction. 
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The beauty of EG?'' process is its ability to machine 
very hard metals without causing and tool wear, high metal 
ranoval rate, no structural damage to the workpiece and to 
machine complex shapes with good surface finish that too in 
single stage. The different fields of applications ' 1 of EC?v! 
are drilling, boring, turning, cavity sinking, milling, grinding, 
deburring, polishing etc. However, the capabilities of such a 
high potential process have not been fully exploited because of 
lack of understanding of the mechanism of metal removal and 
inexistence of an efficient tool design methodology. 

ECl"! is a highly complex phenomenon. Most of the mach- 
ining parameters such as, electric field distribution, electro- 

T 

lyte conductivity (K), gas evolution, temperature rise (<at), 
valency of dissolution (Z), electrolyte flow velocity (u), etc., 
are simultaneously and continuously varying along the electro- 
lyte flowpath and with the machining time. Their interaction 
is of complex nature and there is a lack of clear understanding 
of the mechanism of metal removal. Further, the electrolyte 
in the gap is a mixture of evolved gases (such as hydrogen, 
oxygen and chlorine), solids and precipitates. In majority of 
the cases, lEG is complex shaped and it is difficult to evaluate 
flow parameters since there are electrolyte starvation, and 
cavitational effects as well. 

1.2. TOOL (CATHODE) DESI® FOR ECM 

Tool design in BOA deals with the computation of tool 
shape which under specified machining conditions would produce 
a work having a prescribed shape or profile accurately. 
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Alternatively, it also deals with the production of work shape 
obtainable from a tool of knoivn geometry when machining is 
performed under specified machining conditions. 

There are two main aspects in the design of tools for 
SC?’ _ 2 . First is the determination of the shape of the tool 
together with the optimum machining conditions’_3'^ necessary to 
produce the required work shape. The second aspect is a prac- 
tical one, which is concerned with making the tool of the 
appropriate material, fixing it in the machine, connecting it 
to the power supply, arranging an adequate supply of electrolyte 
between tool and workpiece and insulating parts of the tool or 
to use bit type of tool to prevent machining in undesired 
regions. The design requirements for some of these aspects may 
be conflicting and many of them may make it necessary to modify 
the tool geometry. Here, only the first aspect of tool design 
is dealt with prescribed machining conditions have been used 
for cathode shape design. 

In general, cathode design for ECM is based on the 
‘trial and error’ method which also involves actual machining. 
Such practical methods of obtaining tool shape are expensive, 
time consuming and inaccurate# The analytical design of ECM 
tools is complicated on account of complexity of interaction 
between simultaneously and continuously varying machining para- 
meters. Several attempts have been made to develop a compre- 
hensive model of the mechanics of metal removal in EDA, parti- 
cularly with the view of developing an accurate method for tool 
design* Some of the models for anode shape prediction and 
cathode shape design have been discussed below. 
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1.3. 30S9 -'ETroD 

Tipt©n 4, 5 attempted tool design for EC?1 by using 
cos8 method, with the assumption that the electric flow lines 
are straight and normal to the vioik surface* In this method, 
equilibrium work shape is computed corresponding to the tool 
whose profile has to be approximated by large number of planar 
sections inclined at different angles, as shown in Figure 1.2. 

For plane parallel electrodes with their surfaces 
normal to the direction of feed, the equilibrium gap, Y , is 
given by equation ( 1 . 1 ) 

The equilibrium gap between any section of cathode surface and 

that corresponding surface of the anode which is parallel to 

it is given by Y /cos0. The angle is measured between 

© 

normal to the anode surface and the direction of the cathode 
feed. The cathode shape for a given anode shape can be computed 
from the following parametric equations. 

Suppose the anode shape is given by, 

Y = f(K) (1.2) 

^ny point, say ^(X, Y), on the anode surface corresponds 
to an equivalent point B(X^, Y^ ) on the cathode surface. The 
gap width between A and B then becomes Y^/cos© (Figure 1.2), 
we have, 

Y - Y^ a AB cos© a Yg (1.3) 


and 
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- X = sin0 s= 

By substituting the values of X and Y from equations 
(l,3) and (1*4) in equation (l*2), gives an equation for the 
cathode surface j* 

The *cos9' method is generally applicable to the 
regions where the local radii of curvature of the anode and 
cathode surfaces are large compared with the equilibrium gap. 
This method excludes the effect of different process parameters 
on predicted anode profile. 'Vorkpieces having sharp comers 
and complex shapes cannot be analysed with this method. 

1.4. NUr/BBiaU METHODS 

In E(7% the change in shape of the workpiece during 
machining depends on the local current density, which in turn 
changes with the changing work shape. To determine current 
density, it is necessary to solve Laplace equation for the 
field distribution in the I EG. With the availability of high 
speed computers, numerical methods for solving Laplace equation 
have been developed. 

1.4.1. FINITE DIFFERENCE TECHNIQUE (FDT) 

Tipton [253 snd Hopenfield and Cole [.*6 J applied FDT to 
determine potential distribution in the lEG for non-passivating 
electrolyte assuming constant conductivity and temperature. 

Figure '1.3(a) shown the tool work surfaces and lEG 
drawn in square meshes containing a set of grid points of 


tan© 




(1.4) 
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general co-ordinates (i, j). The potential values for the grid 
points on the cathode boundary are assigned zero, and those on 
the anode boundary are set to some known values. The initial 
potential at the grid points within the IHG region are set by 
linear interpolation along the vertical grid lines between 
anode and cathode. These values must be progressively adjusted 
until they satisfy a finite difference equation corresponding 
to Laplace's equation in the region between the boundaries. 

^or a point 'o' located in a mesh of spacing h{i, j), 
the finite difference equation corresponding to Laplace's 
equation becomes, 


^2 


o 


(1.5) 


Using eq. (1.2), eq. .(1.3),, can be obtained for calcu- 
lating the potential at 'o' 

(0.t , . Jj 0jj . J ■? ■< ) 

J) ^ 


( 1 * 6 ) 


In some cases, all points on th§ anode and cathode boundary may 
not lie on grid points. In such cases, it is Customary to 
relax the potential by linear interpolation. 

1.4.2. FINITE ELB'IHNT METBDD (FEM) 

In case of complex shaped lEG, the cathode and anode 
boundaries cannot foe matched accurately using square meshes, 
which introduce further approximations. In order to overcome 
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these approximations, Jain and Pandeyf?’’ suggested which 
can be applied to analyse any complex shaped lEG. 

Jain' s;^ has developed FH' models, FET11 and FET22, 
for the prediction of anode profiles in the case of plane para~ 
llel machining. The model FET11 is based on the Ohm^s law for 
the calculation of current densities in the I5Q and FH»^! has 
been applied to obtain temperature distribution within the I EG. 

In the model FHT22,' Laplacian equation has been solved to 
obtain potential distribution within the I EG from which current 
densities are calculated. Jain also extended the model FET11 
(Model SGFET11) for predicting the anode profile in case of 
rectangular and cylindrical deep hole drilling operations. All 
these models are capable to predict the anode shape accurately, 
by taking into consideration, the variations in different process 
parameters. ' A new concept of EG?.l tools was introduced by Jain 
[ 9 , kT]. Murugan |jlJ and Yogindra \__i 23 have extended FEhl 
models for anode shape prediction using bit type of tools. 

1.4.3. BOUNDARY ELB'tENT f^ETHDD (BBA) 

Both the methods FDT and FFM yield solutions correspon- 
ding to internal as well as boundary nodes, but in most cases 
(like in EG^A) internal solutions are unnecessary. In BEM, 
boundary alone is discretized into segments. It considerably 
reduces both the amount of input data and its preparation time. 
The solutions are predominantly boundary based but solutions to 
internal points can be obtained. 

Narayananan, Hinduja and Noble [jisj have applied BBA, 
for the anode shape prediction in EQA. Linear and quadratic 
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elements were used to represent the boundaries. It was concluded 
that comparatively better accuracy was obtained by using BSi, 
with far less computational effort. 

1.5. GATHDDE SHAPE PREniGTION BY DRAWING CURVILINEAR 
SQUARES IN THE I EG 

Tipton rsu has developed a method of determining the 
tool shape by dra'.ving curvilinear squares between the work and 
tool boundaries. It defines the ivork surface numerically in 
X, Y plane \<vith Y as feed direction and with coordinates spaced 
by constant X increment equal to a fraction of the equilibrium 
gap. In the equilibrium condition, the normal current density 
is constant, so that a constant current incranent I flows 
between each of the adjacent coordinates on the work boundary 
because they are equally spaced in the X-direction. The worfc 
boundary is assijtmed to be equipotential and each coordinate 
represents the end of a current flow line. By using each adja- 
cent pair of coordinates, it is possible to draw the curvilinear 
squares and to construct the first row of potential field dist- 
ribution between the work and tool as shown in Figure 1.4. If 
the X coordinate spacing is Y /n then repetition of the process 
will produce the nth equipotential which is a tool surface. 

Each of the intermediate surfaces is also a tool surface for 
equilibrium gaps of Y^/n, PY^/n, 3Y^/n, etc. 

The process of drawing curvilinear squares is based on 
differentiating the input data. Differentiation is a noisy 
process, which magnifies random error in the data, specially 
when the process is repeated several times. Complete failure 
of calculation occurs when two current flow lines cross each 
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other. In that case, the previous equipotential could be used 
as tool with the corresponding fractional equilibrium gap. In 
this model, variations in all the process parameters have not 
been taken into account. 

1.6. gOMPLHX VARIAPLr: TECH:JIQUH 

Wilson and Tusei 'H4, isj, used complex variables tech- 
nique to solve the field equation in IE3, The tool design 
problem was formulated as an inverted Cauchy problem, in which 
the spatial coordinates were treated as the dependent variables 
on the plane of the complex potential {21 bj. The equations of 
equipotential lines are directly obtained for a given work 
profile [j^J, which define the equation for a family of tool 
curves. Limitation of the above method is that it can be used 
only for analytic work surfaces, i.e. this method cannot be used 
for discontinuties in derivatives of work and tool surfaces. 

It needs a lot of mathematical analysis. 

1.7. CONTINUITY I’ETIDD 

Lawrence 21^3 ® graphical technique for tool 

design by solving the field equation in the gap. Rut this 
procedure gives low accuracy and it is more tedious. He subse- 
quently adopted numerical methods to solve the field equation 
using digital computer. In this method, starting from a defined 
equipotential work shape boundary, a model of electric field is 
constructed, where each new equipotential surface is a possible 
tool shape which will produce the original work shape. This 
method rests on the assumptions that Ohm*s law applies across a 
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small element of electrolyte and that the current flowing bet- 
ween a pair of flux lines in an electric field is constant. It 
does not account for the variation in electrolyte parameters in 
the IHG. 


1.8. HEAT TRANSFER IN FB! PROCESS 

During ECH, a large amount of heat is generated in I EG 
by different sources of heat viz.. Ohmic resistance of electro- 
lyte, reversible and irreversible chemical reactions, viscous 
flow of electrolyte etc. A major contribution is by the Ohmic 
heating of the electrolyte |17]. But a significant amount of 
heat is transferred to the surrounding atmosphere through the 
electrodes. This lowers down the temperature of the electrolyte, 
hence it becomes essential to incorporate the effect to modify 
the electrolyte conductivity along the flow path to get accurate 
results. Apart from the conductivity, change in electrolyte 
temperature affects the other machining parameters also. Heat 
transfer analysis in EG*' is quite complicated and involves 
multiphase turbulent flow with heat generation between the 
electrodes. In majority of the models for ECM tool design, 
therefore, heat transfer through the electrodes has been negle- 
cted. This leads to poor correlation between analytical and 
experimental results. 

1.8*1. THERf.lAL RESISTANCE ,TODEL 

Jain and Pandey have suggested a simplified, 

unidimensional thermal resistance model for heat transfer anal- 
ysis of ECfi^. For an IHQ system, temperature gradient and 
simplified thermal resistance model \^18'J are shown in Figures 
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1.5(a) and 1.5(b) respectively. The drop in temperature of the 
electrolyte due to transfer of heat to the atmosphere is eval- 
uated using the folloiving equation, 

q *= Qa + 

where q_ = heat transferred from anode to the air 

and q = heat transferred from cathode to the air. 

qut this method is based on the following assumptions; 

1. IHCl system behaves like a single phase fluid 

2. Heat transfer through electrodes is unidirectional 

3. Heat transfer within the occurs under steady state 
conditions. 

Murugan Q11, I9j had investigated into the effects of 
cathode material on temperature distribution during EC?';, using 
the simplified thermal resistance model. 

Kawafune, Mikoshiba and Note analyzed the effect 

of tenperature and other working parameters on the working gap 
during EQ/.. In the analysis, an equation for the temperature 
rise had been derived assuming steady state heat transfer. Also, 
the physical constants (K, P^, ) of water were used for 

heat transfer equations instead of those for electrolyte* 

1.9. PRESENT WORK. 

From the above literature survey it is evident that 
there is no accurate model available for heat transfer analysis 
of EC*''. The work embodied in the first part of this thesis has 
been aimed towards the development of a two dimensional, unsteady 
state heat transfer model to analyze the heat transfer fromth 
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the electrolyte to and through the electrodes, for obtaining 
the temperature distribution in the I EG. The anode shape pred~ 
iction model (FET22) for plane parallel electrodes HCM developed 
by Jain i"8' , has been modified by incorporating the two-dimen- 
sional, unsteady state heat transfer. It has been named as 
TFET22. 

The analytic temperature distribution, obtained by 
TFET22, in lEG has been compared with the experimental tempera- 
ture distribution obtained by Bhatia [ 21, 22j and with one- 
dimensional heat transfer analytical results obtained by 
Wurugan (”*1 1 ' ( . 

The second part of the work has been aimed towards the 
development of a computer aided method for the tool design for 
ECf/. The cathode (tool) shape prediction problem has been 
attempted for electrochemical drilling with a bare tool, electro- 
chemical bit drilling and electrochemical bit boring. The anode 
shape prediction model (SBFET22), first developed by Jain (^8 j 
and later modified by Yogindra '[12’j, has been used in the design 
package for ECO and ECBD. The model SBFET22 has bem modified 
for the anode shape prediction for ECB. 

Tools and tool bits have been designed for producing 
profiles obtained experimentally and also for producing straight 
sided holes, while machining under prescribed machining condi- 
tion's. hesigned tool shapes have been compared with the experi- 
mental tool shapes used by Jain for ECD and ECBD. 




FIG. 1.1 BASIC SCHEME OF ELECTRO-CHEMICAL MACHI- 
■ NING PROCESS, 




(a) Principle of the cosB method (b) Cathode design by cos 0 method* 


FIG. 1.2 Cose METHOD 
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CHAPTER II 


THEORETIC \L AH-ALYSIS 


2.1. IMTRODUCTIOH 

In this chapter, EC.Y with plane parallel electrodes 
and electrochemical drilling have been analyzed separately. 

The expressions for various process parameters such as electro- 
lyte conductivity, flow velocity, current density, void frac- 
tion, temperature and pressure distributions etc. have been 
given. All the three zones namely, front, transition and side 
zones in the case of "Ul have been analyzed individually with 
some suitable assumptions. However, stagnation zone has been 
igno red. 

2.2. ANAI.YSIS OF F.a^! ”/ITH PLAI'IE PARAIXEL ELECTRODES 

2.2.1. Ef , ECTROLYT E FI..0'7 VEIOCITY 

Assuming uniform velocity distribution across the flow 
direction, electrolyte flow gielocity at any cross sectiton is 
given by 

U *= I (2.1) 

where Q is the discharge and A is the cross sectional area. 

2.2.2. CURRENT DENSITY AND CONDUCTIVITY 

The current density within the I EG, can be determined 
either from Ohm's law or from the potential distribution obtained 
by solving Laplace equation. The application of a particular 
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case depends upon the nature of the current lines. Ohm's law 
can be applied only when the current lines are straight and 
parallel. But this is not true when electrodes with small 
radii or complex shape are used. In such cases, the current 
density calculated from the potential distribution obtained by 
solving Laplace equation gives more accurate results. The 
current density at any point on the workpiece is calculated as 

J = K (2.2) 

(,.n 

The electrolyte conductivity, K, has been known to be 
a function of both, electrolyte temperature and void fraction. 
For computation of K, Tipton proposed the equation given 
below, which accounts for the effects of electrolyte tempera- 
ture and hydrogen gas liberation: 

K = K :i +<At‘|( 1 - (2.2a) 

In eq. (2.2), Hopenfeld and Cole []6j assumed a value of n = 1.5 
for the case of uniform void distribution. For the case of 
non-uniform void distribution, especially where bubbles are 
concentrated near the cathode, Thrope and Zerkle ^23j suggested 
a value of n « 2.0. 

2.2.3. CURRENT EFFICIENCY 

In deriving the expression for metal removal rate (MRR) 
by Faraday's laws of electrolysis, current efficiency was 
assumed to be 100®^. But in the actual machining process, it 
depends on the current efficiency achieved (^l]]t i»e., the 
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proportion of the total current that is used for the removal of 
metal from the anode. The current efficiency is given by 

'A - Actual metal removed (o 

• ” Theoretical metal removed ’ v • ; 

In actual practice, the values of current efficiency 
may be higher or lower than 100'^'^, depending on the conditions 
of machining. Current efficiency values less than 100’'' are 
obtained because of the occurrence of the side reactions at 
the anode, e.g., evolution of oxygen gas or Fe ions instead 
of Fe'^**’ ions etc. 

If the dissolution of the anode occurs at a valency 
lower than the one used during calculation, the current effi- 
ciency obtained may be higher than 100'^’. Voreover, selective 
machining may occur at the material irregularities, such as at 
the grain boundaries and hence may cause chunks of material to 
be r©T50ved. In this case, the actual material removed is more 
than the theoretical which results in current efficiency higher 
than 100^<. 


2.2.4. IMTFRFLFCTRODH W (IHG) 

TG reactions would attain an equilibrium only after 
time t approaches infinity and IFG can be computed from equations 
(2.4) and (2.5) while machining with zero and finite feed rates 
[^1*j respectively, 

Y *= (Yq 4 2ct)‘*'^^ (2.4) 

. AK(V - AV) 

where c » 


and 
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y — y ~ 

t = (y^ - y) + ' (2.5) 

— e - 

Equation (2.5) is an implicit equation and can be solved by 
iteration process. Equation (2.4) predicts an infinite value 
of the gap as ' t’ approaches infinity. However, in practice, 
as the lE'l increases, the current density decreases, so that 
the rPJl gradually decreases and for a very large gap, the 
process would come to a stand still. 

Eor a very small interval of time !it, Jain and Jandey 
1^7, 10J derived an oq. (2.5) for determining lii which holds 
good for both, zero and finite feed rate. 

Y = y^ + (c* - f)At (2.6) 

V- , J A 

where c’ = •■• y - p- . y 

Equation (2.6) is based on the assumption that the 
current density ’J' remains constant over a small element of 
length dX and for a small interval of time, ^t. This equation 
yields values of Y which are close to those obtained from eqs. 
(2.4) and (2.5) and agree well with experimental data[[8'J. In 
the present work eq. (2.6) has been used. 


2,2.5. TEYPHXATUXE 

Following assumptions are made in deriving the equations 
for temperature distribution! 

(1) Total heat generated is only due to Ohmic heating (i.e. 

I'^R^), The heat generated by other sources is negligible l‘l. 

(2) Specific heat and density of the electrolyte remain cons- 
tant in I EG. 



20 


(3) Heat transfer through the electrodes is considered to be 
negligible i.e., all the heat generated is retained by 
the electrolyte. 

Consider an element of length dv and width in the 
IBG (Figure 2.1 ). 

The mass of the electrolyte flowing 



2.3. WAiYSIS OF FC ORILUNG 

2.3.1. FLHCTROLYTE FLO'V VELOCITY 

In the present case, the mode of electrolyte flow is 
outward and separate expressions have to be derived for front, 
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side and transition zones* The assumptions made are, 

(1) The volime of electrolyte flowing at any section per unit 
time is constant. 

(2) Cavitation and starvation of the electrolyte are neglected. 

(3) Electrolyte velocity across the flow direction is assumed 
to be uniform. 


7710MT ZOr'IE 

Tn this zone, electrolyte flows in radial direction. 
Electrolyte flow velocity at any radius (Figure 2.2) from the 
axis of the tool is given by 


^ 2TVrY 

where Q « discharge of electrolyte. 


(2.11 ) 


TMNSITIO^^ AND SIDS ZONES 

As the tool comer radius is very small electrolyte 
flow in both the zones is assumed to be in axial direction. 
From Figure 2.2 

= TTCr^ + Y)^ ^ Ifr^ 

= Tr(Y^ + 2r2 Y) (2.12) 


and 

U * "g" - 

K (Y^ + 2r2 Y) 

2.3.2. CiJRRENT DENSITY AND CONDUCTIVITY 


(2.13) 


These are evaluated in the same way as discussed earlier 
in the case of ECM with plane parallel electrodes. 
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2.2.3. TS'?3R\PRE 

By making the same assumptions as in the case of plane 
parallel electrodes the expressions for temperature distri- 

butions have been obtained as follows. 


PROMT ZONE 


Consider an element of thickness *dr’ at a distance, r, 
from the axis of the tool (Figure 2.2). 

Using Ohm* s law, it can be shown that the current I 
through this element is given by 


I 


K (V - Ay) 2'V. r dr 


(2.14) 


Heat generated in the element. 


dl! = I^R. 


:(V - AV)^ 2 T r dr 
^ 


(2.15) 


Heat carried away by the electrolyte, 

dH = 0 dT 

e e 

Equating (2.15) and (2.16) yields 

Q p c dT = (V - AV)^ K 201 r dr 

e e ' '' ^ ' 


(2.16) 


(2.17) 


Using K = Kq( 1 +a<AT), the expression for temperature distri 
butidn becomes. 


^ (2.18) 

♦ 2 K„(V -AV)^ 

where, A =:= — 

* e e 
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The integration of eq. (2.18) with boundary condition at 


T = T. gives. 


T - Ti 


exp-J A 


f 2 2 \. 

(r - r^) 


(2.19) 


TRANSITIOM AND SIDH ZaM:-3 


Jn these zones, the electrolyte flow is assumed to be 
in axial direction. Consider an element of length dX at a 
distance X from the end of the front zone, i.e., from 
(Figure 2.2). 


Oy Ohm' s law, 


(F -.AV)K 21l(ro + Y)d'' 


■rom eq. (2.20) the heat generated in the element is given by 


(V - V)^K 21l(r2 + V)dX 


( 2.21 ) 


■quating eqs. (2.16) and (2.21) gives 


0 2^ dT 

e e 


(V - aV) K 2iT(ro + Y)dx 


( 2 . 22 ) 


A’ dX 


where. 


2Ti K (V - AV)^ (r^ + Y) 

cTfTv^ 


(2.23) 


Integrating eq. (2.23) with boundary condition, i.e. 


^ref “ ’'ref 


(T - T, ) 


exp|'< A’(X - 1 : (2.24) 
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where, 

'3 


1 +c, (T 


rei 




(2.25) 


2.3.4. VOID FR.^CTION 

During high current electrolysis, hydrogen will evolve 
at anode and oxygen the cathode. Hopenfeld and Cole are 
of the opinion that the amount of oxygen is too small for real- 
istic considerations. Similarly, the presence of chlorine can 
be neglected due to its solubility in water. Hence, for void 
fraction calculation, the effect of hydrogen is considered. 

Considering a control volume in IHG (Figure 2.3), 
using the principle of conservation of mass, Thrope and Zerkle 
[23' arrived at the governing nonlinear differential equation 

* „f0* 

and the expression for void fraction was obtained by solving it, 
as follows 


V 


■!f * 

B X 


1 -t 3 


(2.26) 


* 

where, B 


♦ 

X 


s 


& 

' e 
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(for rectilinear flow) 


X(r. + X/2) 

— (for radial flow) 

L o Ti 


Pi "i 

TTrir 


In the present case, the gas slip ratio <r' in the variable 
B* is assumed to be unity, that is, velocity of gas is same as 
that of velocity of the electrolyte. 
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2.3.5. HT.^CTHOLYT^ PR^SSim 

The electrolyte pressure "1^ required to maintain the 
rate of electrolyte flow in a particular gap is given by 


p, + 


(2.27) 


where is the pressure required to overcome inertia force 
and p 2 is the pressure required to overcome viscous force. The 
expressions for and ^2 different for laminar and turbulent 
flows. The type of flow in the 133 can be determined from 
Reynold's number and is given by the eq. (2.!S) 


U D 
© © 


(2.28) 


where D = The hydraulic mean diameter 

_ 4 X Area of cross section 

“* 'Vetted parameter 

= d for a cylindrical tube 

“ “ *^2^ case in which fluid is flowing 

through the annulus 

s= 2Y for a rectangular channel of depth Y. 

the laminar flow (^41 2000) with a parabolic velo- 
city distribution, the pressures and P 2 are given by 

P. = i f 

1 2 ' e e 


32 U X 

■ " ■■■ "K - — for a narrow tube of 

d'^ diameter d 


(2.29) 


The flow is usually turbulent (3^> 2000) in E® 
From Bemouli’s equation, 
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P 


1 



(2.30) 


The value of is given by Dlasius ecuation 

0.3164 i'\ 


2 

f *f 

_ Vi _ .V 

e e 


n 


The pressure Pr,, i?^ turbulent flow is a function of 
Reynold’s number. At hiqh Reynold's number inertia force will 
be more than the viscous force which results in higher value of • 


2.3.6. T-'l 

TRl is calculntod in the same way as discussed earlier 
in the case of IT' ^vith plane parallel electrodes. 


2.3.7. FERH RATE IN TRANSITION AMR SINH ZONES 

In EC drilling, as the machining progresses, the side 
zone becomes tapered and it no longer remains zero feed rate 
case, ^or the calculation of IHG in the side zone, the feed 
rate vector normal to the work surface should be considered, 
as given by eq. (2.31) 

fj^ = f cos«5 (2.31) 

where is the angle of inclination between the feed direction 
and normal to the work surface. The values of ’9’ for each 
element in side zone and transition zone are calculated separa- 
tely. 

3f Y. and are the interelectrode gaps, say in side 

1 Z 

zone at points 1 and 2 resnectively (Figure 2.4) then the angle 
9 at the point 2 is given by 
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1 

^ = tan”^ (7^ y^) 

In the transition zone, the anqle 9 at any point is 
given by the intersection angle between the normal drawn at a 
point and the tool feed direction (Figure 2.4). 



2S 


r 


ELECTROLYTE 



FIG. 2.1 PLANE. PARALLEL ELECTRODES ECM. 



ELECTROtyfC now AREA 



(a) CURRENT FLOW AREA 


SIDE ZONE 





Id) FRONT ZONE (a) TRANSITION AND 

SIDE ZONE 


pig.2.2 Electrolyte and current flow area for outward mode 
of electrolyte flow. 
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CHAPTER III 

HEAT TRANSFER ANALYSIS OF ECr PROCESS - f.^ODEL TFET22 

3.1. INTRODUCTION 

The aim of present day research in EC?' is to devise 
an accurate and realistic methodology for tool-design. Conven- 
tional ECF^' tool-design is not accurate enough to account for 
heat transfer through the electrodes. Due to the complex 
nature of the process, in majority of the models for ECU tool- 
design, heat transferred from electrolyte to the electrodes has 
been neglected. 

Tool materials generally employed during EOU are copper, 
brass, stainless steel, etc., which are good conductors of heat. 
Since large amount of heat is generated due to the Ohmic resis- 
tance of the electrolyte, considerable amount of heat will be 
transferred to the surroundings through electrodes. As the 
conductivity of the electrolyte, which is the most important 
machining parameter during ECM, is very sensitive to the temp- 
erature. Therefore heat transfer during EC^' demands in-depth 
analysis. 

Jain and Pandey have suggested a simplified thermal 

resistance model for heat transfer analysis during EO^i based on 
simplified assumptions like, unidirectional heat transfer, 
steady state heat transfer etc. It was concluded that the 
effect of heat transfer on the machining parameters during ECM 
is significant and ECM process must be analysed as a two or 
three dimensional heat transfer problem. 
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Here, an attempt has been made to analyse the heat 
transfer during HC” as two dimensional. Unsteady state, con- 
jugate heat transfer problem. Finite element method has been 
applied for solving the governing differential equation. 

3.2. CONJUGATE HEAT TRANSFER 

When a solid material of finite dimensions and mode- 
rate conductivity is adjacent to a fluid flow domain, the temp- 
erature distribution in the fluid is coupled with the temperature 
field in the solid. Often, the thermal boundary condition is 
not prescribed at fluid-solid interface, but it is prescribed 
on the outer surface of the solid. This represents a conjugate 
heat transfer problem. If the conduction in the solid and the 
convection in the fluid were to be solved separately, then the 
matching procedure at the interface would have been very diffi- 
cult. Here, the finite element method offers an easier altern- 
ative. 

Since the finite element method can admit discontinuous 
conductivity distribution, all that is needed is to specify the 
solid conductivity to all the elements in solid region and 
fluid conductivity for all the elements in the fluid region. 

The solid-fluid interface will coincide with the control volume 
face, hence no boundary condition is needed there. The specified 
boundary condition at the outer surface of the solid can be 
directly used at the boundary the computational domain. The 
resulting solution will include the temperature distribution in 
the solid as well as the fluid. 
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In EC?' process the electrolyte flows between tool and 
work (Figure 3.1). The heat transfer problem during ECM is a 
perfect conjugate heat transfer problem. The finite element 
discretization of the computational domain is shown in Figure 
3.3. Once the governing differential equation is derived, the 
temperature distribution can be obtained by the finite element 
method. 


3.3. FOm^ULATION OF THE PROBLB.' 


To derive the governing general differential equation, 
consider Figure 3.2 where energy balance analysis has been done 
for an elemental volume in the machining system as follows s 


Energy in + Energy generated = Energy out 
- K 4 Y - K AX + f Cp U T A Y + H 




+ f C M(U +'^AX)(T +f^Ax)4Y 


Simplification of about equation leads to 



ss 



f 




H 


(3.1) 


The above equation is derived assiming unidirectional 
fluid flow (x-direction, only) and negligible heat generation due 
to viscous flow. To get the temperature distribution ^ it has to 
be solved subject to the boundary conditions as in Figure 3.3. 
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Equation (3.1) is fully applicable for an element in 
liquid zone. For an element in the solid region, the convection 
and heat generation teims will be zero. 


3.4. FIMITE ELB-’ENT ANALYSIS 


The governing differential equation is, 

2„ 


K 


^ . '"T 

I K H- » 

» X , 


^ P r ij “iX 4 . 


H 


r- '"Lx 

' "p f,: t 


0 


(3.2) 


’(Vith boundary conditions 

(a) Temperature may be specified along the boundary 

(b) It may be a convective boundary. 

Let the temperature variation in the element be as follows: 


T^ s= a + bx + cy 



Then, the residue (R) is 


K 








^T , 

P C U 4 - H 

'' P ^x 


K *r* 

Ifi i 

A n 


(3.3) 


Gale^dcin method is used to minimize this residue. 


. • . ) Nji n dx dy « 0 

(3.4) 

First term of the eq. (3.4) is 
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, 2^ 

K dx dy 


Integrating with respect to x *24 gives 


r; ■! , 


i<. N. — - A- dx dy 

t X ♦. ^z. 

tJ X 


K N 


e 


1,X d X 


dx dy 


+ ■ K N. D dB 

. i,B X X 


where is direction cosine .along x-axis on the boundary of 
the domain. 

Integrating the second teim with respect to y gives 




! 1 > 2 
^ uy 


^ dx dy 


^i,y if 
^ ^i ,3 zf 


where n^ is direction cosine along y— axis on the boundary of 
the domain. 






r)T 


(( kI’^ + ^'1 dx dy = - (< KfN 'r-^ + M r-^ 


^ y 


jdx dy 


<^T, 


JT, 


+ ' ►'ibC*^'^ + K r-rrS n.. I dB 


'1B\_ jx “X ’ "jy y 


"yl 


fi + My 5^'Jdx dy - Jp 

fi'K [{n;^ LNxi-^l^yJ LNyj]dx dy(Tj- J h{T - dB 

- W{'f] - ” U'^Bis l-Nj fl^ dB + h T„ f Nb dB 
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Third term is. 




0 u ' ■ .'n. ' N ' dx dy/T , 


A ' . t'. 


Fourth term is, 

H dx dy 

- H 5 dx dy = P 

Fifth term is, 


'p / ‘ dx dy 


I'S i'5\N} UJ dx dy^i;. 


.*» ‘ 1 » M 

-! M 1 1 T '1 


Finally by combining, 



-Ik' 

w 

: ;;'t‘'. - ' Bi' ^ T ♦ -f !'B2 i 

* \ ■/ J w 

or 


l®il ■*■ [^.3 + I'-i 

Let 

[B] 

= [Vj - [®0 ^ 'A1 

and 

I'BJ 

= {pj + ^B2'3 

Then 




+ N P t - '! M ' % T 






<3.6) 
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3.4.1, GHMER^TION OF DIFFERENT MATRICES FOR TRIANGULAR ELBIENTS 

For triangular elements (Figure 3.4(a)), the t®itperature 
is assumed to vary as 

T = a + bx + cy - 1 . 

Let 



Area of the element is given by, 

A .= i [ x,(y2 - Y 3 ) + - yp + "a^yi - y20 

Interpolating functions are given by, 

= k - ’'2y ^ 2^3 - h '^2 + - ’'ysj 

*^2 = k 1*1 y - *yi + *y3 - *3^ *3yi - ^lysj 
^3 ' k 1*11^2 - ='2i'i + *2y - *y2 + *yi - ’'ly'] 

“ k 1*^1 ^2 

% 

^,y 2A l?1 SJ 

The matrix [kJJ is 

t^s) = U [('^xU'^xi +[\yLvJ'’* 
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4A 


, 2 2 

+ of 

1 1 


k2 . 2 

^2 <^2 


Sym 


' A 


«? O T T 


1 N * . N _ dx dy 

*■ 5v 


'ITi r» |T 

> V 


lif 


b. 


I P 


, ' M ■■♦ dx dy 


HA\ 


f : 1 

V 1 j 


^ 1^3 ^ 1^3 

^ 2^3 ^ 2^3 

k2 ^ ^2 

^3 ^3 


3.3.2. BOUNDARY MATRICES 


The boundary matrices are: 



[ B2 


h ( ( N > ! N • dB 

B 

h T^ \ N% dB 
o B ^ ^ 


Thus for evaluating these convective boundary matrices, 
the convective boundary B must be known. If the edge i,j is 
along the convective boundary (Figure 3.4(a)), then. 
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li • • S • » 


If the edge jk of the element is along the convective boundary. 


then 


h»S ji 


0 0 0 
0 2 1 


0 1 2 


h.T_.S, 


If ld[ is along the boundary, then 


. h.S,, 


h.T .S., \ 

Cp Kk X 0 '*'1 

I ' . 


3.5, METHOD USED FOR SOLVING THE SYSTB'I OF LINEAR DIFFERENTIAL 
EQUATIONS 

In order to get the temperature distribution, the system 
of linear differential equations given by eq* (3*6) has to be 
solved. The following method is adopted for that purpose. 


.{b'3 


(3.6) 


From the Figure 3.4(b) T is given by 
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* k+l > T « 

> T« - 9 i' t (Backward integration) 




f T, + (1 - 9) rrr -■'■‘t (Forward integration) 


where 


© is a fraction 0 9 1 


t is time interval 


V; ' - ©(T)^^'"'^^. t 


M ' I , 4(1- 9)/T'.^ At • 

' i, • 

(3.7) 


From eg. (3.6) 


^T* \ k4 1 


r 4 s n' -rrr^' * v b' 


■‘T-'i J«4l 
'1^1 


' M A' T 'A 4 0 . “ T * 


Using egs. (3.8) and (3.9) in eg. (3.6) 


(3.8) 


(3.9) 


[(^j n + S At “J 

« 4 (1 - 9) Atr^B^^ - [dJA* Tj^J 


For ) B f 


1^/ 


+ 9 At [D]] = ^[M] - (1 - 9) At i + 

[^(1 - 9) At{ B'j 4 9 ^t{B^ '3 
For 9 a= 1 (Backward integration) 


Ml + AtfollfTA"-' 




V"j{U +At(BJ 


Dividing by A t 


AW-i'U 


■AeD-0(Tf 


(3.10) 
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Now the problem is simplified. To get the temperature 
distribution, the system of linear equations given by eq. (3.10) 
has to be solved. 

3.6. RESULTS \ND DISCUSSION 

Electrochemical machining experiments with plane parallel 
electrodes were conducted by Bhatia ^ 21 , 22 using copper, brass 
and stainless steel as tool materials, mild steel as work 
material and NaCl solution in water as electrolyte. Temperature 
at distances as 2.5, S.fj, 17.5 and 21.5 mm from one end were 
measured using thermocouples. The experimental conditions and 
temperature data are given in Table 3.* Uurugan [^1 lj conducted 
one dimensional heat transfer analysis of ECM process using 
thermal resistance model tested his model using experi- 
mental data reported in reference this analysis, heat 

transferred by conduction mode alone was considered, and heat 
transfer by convection was neglected. As a result of such 
simplified assumptions a discrepancy (Figure 3.6) is observed 
between analytical and experimental results. Keeping this in 
view, the problem has been analysed as a two dimensional, single 
phase transient heat transfer problem. 

Figure 3.6 shows the comparison of experimental results, 
one dimensional heat transfer results and the results obtained 
by the present model for different tool materf.als. Results have 
also been plotted (Figure 3.5) for the case, when heat transfer 
through the electrodes is neglected (Model FET22). It can be 
seen from the figures that the results obtained by the present 
model are closer to the experimental data. This clearly shows 
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the importance of incorporating two dimensional, transient 
heat transfer in Hr;?' analysis. However, the discrepancy between 
experimental and analytical results can be further reduced by 
analysing the problem as a three dimensional multiphase tran- 
sient heat transfer problem. 

Figure 3.5 shows the temperature distribution along 
both the tool and workpiece. It can be seen that, for stainless 
steel material, the ternoeratures along the tool and workpiece 
are almost equal. But for copper as tool material some differ- 
ence is obseirved between the temperatures along the work and 
tool. This is mainly because of different thermal conductivi- 
ties of work and tool materials. 

However, comparatively large discrepancies between the 
analytical and experimental results (Figure 3.6(d)) are exhibited 
at the inlet side for the case of copper as tool material. 
However, the deviation in case of present model is less when 
compared to one dimensional heat transfer analysis. This is 
because of the fact that there will be some transfer of heat 
from exit to inlet side due to temperature gradient. 

The discrepancy between analytical and experimental 
results can be attributed to the following reasons. Effect of 
void fraction has not beei considered in this model. Exact 
values of valency of dissolution and machining efficiency are 
not known. The sources of heat, other than Ohmic heating have 
also been neglected. 

3.7. CONCLUSIONS 

From results and discussion, following conclusions have 


been drawn* 
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(l ) The two dimensional, transient heat transfer model yields 
results that are closer to the experimental data when 
compared to the results obtained by thermal resistance 
model or one dimensional FH model. 

(2) Cathode material has significant effect on temperature 
distribution along the electrolyte flow path. Consequ- 
ently it would have an effect on machined profile of the 
component. 
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FIG. 3.1 SCHEMATIC DIAGRAM OF ECM PROCESS. 
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CHAPTER IV 
TOOL RESIGN 

4.1. IHTltODUCTION 

The most important aspect in the design of ECV tools 
is the prediction of tool shape to produce the specified work- 
piece shape* It is well know that, tool shape is not simply 
the workpiece shape reduced (or increased) by some machining 
gap dimensions. The gap is not constant over the machining 
surface. It is a function of different machining parameters. 

Cue to complex nature of ECH process, majority of the 
models proposed for EC*' tool design, do not account for the vari- 
ations in the process parameters such as t®nperature, void frac- 
tion, etc. For an accurate tool design in EG?’, variations in 
all the process parameters and their effects on anode profile 
should be considered. 

’Correction Factor Method’ has been suggested by Reddy 
26 for the purpose of tool design in EG?'. In this chapter, a 
two dimensional model for predicting a cathode configuration 
which c3n produce specified anode shape under prescribed machining 
conditions has been described. Cathode shape design problem has 
been attempted for different cases viz., ECD, ECBD and ECB. 

Anode shape prediction model, developed by Jain and Yogindra 1Q12 
has been used for ECD and ECBD. This model has been modified 
for the prediction of work shape with a tapered tool. This 
model has been used in the tool design' package. It has also been 
modified for the prediction of anode shape in ECB. 
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4.2. GFJ>!’~RAL 2ESI3N PuOOEnURE 

Tool design procedure for producing the required work 

shape while machining under given conditions is summarized as 
follows# 

INITIAL TOOL SMAPE 

The first step is to assume a tool shape depending on 
the required work shape and type of 33' ’ operation e.g. ECD, 

ECn etc* For most of the cases, the tool shape can be 
assigned as congruent to the desired work shape. 

’■.VORK SHAPE PREDICTION 

The second step is to predict the work shape obtainable 
from the assumed tool shape* There are many models available 
for anode shape prediction in ECM. The models based on numerical 
methods like FE?1, can predict work profiles within a reasonable 
accuracy. 

ERROR AND CORRECTION 

The third step is to compare the predicted work shape 
with the required work shape and to obtain deviations betweai 
them. If the deviation is more than the specified tolerance, 
correction required to minimize this error has to be calculated. 

NOniFIOATION 

The fourth step is to modify the present tool shape by 
applying the correction appropriately. Using this modified tool, 
again the anode shape is predicted. Steps 2, 3 and 4 are 
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repeated till the work shape within the desired tolerance is 
obtained (see Figure 4*4(d)3, This procedure consists of a 
number of design cycles* Tool is refined progressively in each 
design cycle. Finally it gives the tool shape wlithh will pro- 
duce desired work shape. 

From the above procedure, it is clearly evident that 
the soul of this method is anode shape prediction model. f\ny 
inaccuracy associated with anode shape prediction model will be 
reflected on the designed tool shape. 

4.3. HLECTFiOCHB^.ICAL DRILLING 

Here, an attempt has been made to design, 

(1) Bare tool for electrochemical drilling 

(2) Tool bit for electrochanical bit drilling (EGBD). 

4.3.1. INITIAL TOOL SHAPE 

Initial tool shape has been assumed as cylindrical one, 
having a radius equal to 1 mm lesser than the required hole 
radius at the end of transition zone. In case of transition 
zone, tool comer radius has been assumed again as 1 mm lesser 
than the required work comer radius. 

4.3.2. PREDICTION OF ANODE SHAPE 

In the anode shape prediction model, the field distri- 
bution is obtained by solving the Laplace equation (-~4 *= 0) 

«x ay 

with the following boundary conditions, 

(a) w 0 along cathode (tool) 

(b) 0 « V - AV along anode (woik) 

(c) » 0 along open edges. 
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Since is a moving boundary problem, the domain 
length and yd-dth will change after each computational machining 
cycle time, ,',t sec. As the tool penetrates in the work surface, 
the total number of nodes would also increase# The coordinates 
of each node, therefore should be adjusted. Figure 4.2 shows 
renumbering of nodes after one cycle of computational machining. 

In '< t sec., tool penetrates to a depth equal to the length of an 
element which is assumed to be same in all zones. Here, an assum- 
ption has been made that the feed to the tool is given only 
after machining has taken place for a time interval of/^t sec., 
i.e., the process is considered like intermittent cutting rather 
than continuous one, as in actual HC?.’. 

The IFG’s in the beifinning of a computational machining 
cycle will get modified in the following cycle. Also, the tool 
will occupy new position after a computation machining cycle 
(Figure 4.3). The modified lEG’s are assigned to the correspon- 
ding points on the tool in new position, in order to obtain work 
shape at the end of that computational machining cycle (Figure 4*3). 

4.3.3. CALCULATION OF ERROR AND CORRECTION 

The coordinate system for all the three zones viz. , 
front, transition and side zones, is shown in Figure 4.1. The 
coordinate system is selected in such a way that the error 
calculation principle remains the same for all the three zones. 

The error and correction at each node are calculated as follows: 

(4.1 ) 


ERRORd) » O'.VS(I) - RWS(I) 
CORR(I) « - ERROR(I) 


( 4 • 2 ) 
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where 




o’lsd) 


h-';s(T) 


= 'Correction to be made at Ith node 
s= Y coordinate of the obtained work shape at 
Ith node 

s= Y coordinate of the required work shape at 
Ith node. 


The correction need not always be equal to error in all 
the design cycles. It can be calculated from the trend of the 
past design cycle. This nay some times give less CPU time. 

7rom the second design cycle onwards, the correction can be 


calculated as given below: 




COHRk^ld) X iiRRORkd) 

- ERROi^rrj 


(4.3) 


where suffix k stands for the present design cycle and suffix k-1 

stands for the previous design cycle. 

Equation (4.3) gives the correction that has to be 
applied (TOUI^(I)) for minimizing the present error (ERROFI^Cl)). 
In the present wo3± eq* (4.2) has beei used. 


4.3.4. ’'.O^IFICMIOM 

The tool used in the previous design cycle is modified 
by applying correction to the Y coordinates of the tool at 
different nodes. However, in the side zone, the correction 
should be applied in normal direction to the work profile at 
that node, since the machining takes place in that direction 
(Figure A.2. Appendi:<-l). But in the anode shape prediction 
model, lEG’s are considered to be perpendicular to the axis of 
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the tool. correction is also applied along the same direction, 
in order to obtain new tool shape. 

NIS(I) = Y(I) + CORRd) (4.4) 

where MTS(I) » Y coordinate of the new tool shape at Ith node 

yCI) =: y coordinate of tool in previous cycle. 

Correction at any node (CoriR(l)) may be positive or negative 
depending on the deviation betv;een the ontained work shape and 
required work shape at that node. 

In HCD, unlike plane parallel electrodes EOT, as the 
machining progresses, a node on workpiece is exposed to different 
nodes on the tool at different times. This makes designing more 
difficult in HCD with bare tool. If the tool shape at a node is 

changed, that it will affect the work shape at many nodes, on 

its way. Later problem is eliminated, to some extent, in case 
of HCBD. 


4.4. TOOL BIT DRSIGM FOR HCBD 

Beyond the bit height region, sides of the drilled hole 
or predicted hole, during HCBD, will be reasonably straight. In 
the bit height region, it will be tapered. Tool bits can be 
designed to produce straight sided through holes with a reason- 
ably good accuracy. Initial bit shape is assumed to be 1 mm 
smaller in diameter than the required work shape. The designing 
procedure is on similar lines as discussed above. Only differ- 
ence is in the anode shape prediction model, that is the nodes 
beyond the bare bit height will remain inactive. 
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4,5. HLECTR0CHEtl3\L BIT BORING (EGBB) 

In borinq, the aim is to aniarge the predrilled hole 
and at the same time straighten its sides (or bring vdth uniform 
diameter along its axis) and finish it by correcting any other 
error in the profile. 

4.5.1. anode shape PRHDlrmON 

During Ecre, the following three types of predrilled 
holes can bo encotjintered (Figure 4.4): 

(1) No interference (or clearance type), i.e. predrilled hole 
diameter is bigger than the bit diaraeter throughout its 
depth (^inure 4.4(a)). 

(?) Full interference, i.e* predrilled hole diameter is smaller 
than the bit diameter throughout its depth (Figure 4.4(b}). 
(3) Partial interference, i.e. predrilled hole will be bigger 
upto certain depth and for the remaining depth the hole 
diameter will be smaller than the tool bit diameter 
(riqure 4.4(c)}. 

To calculate the initial I EG for different nodes, the 
predrilled hole profile has been approximated by a polynomial 
curve. Since the sharp cornered tools were used in EGBB process, 
the transition zone during the computation of anode profile of a 
bored hole has hem ignored for the first type of predrilled 
hole (clearance type). Analysis becomes quite complicated for 
the case of interference type of predrllled hole, involving all 
the three zones* For those cases, transition zone has been 
considered from the interference point onwards. Front zone is 
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treated in the same way as in ECD. The other details have been 
given in 4ppendix-I. 

4.5.2. T'lOL HIT DESIGN FOR ECTB 

T'ool hits have been designed to enlarge and correct the 
predrilled hole. The initial tool bit is assumed to be 1 mm 
smaller than the desired hole. The designing procedure is on 
same lines as in ECBD. The other details have been discussed 
in ‘\ppendix-I. 

4.6. RI-Gl.'T.TJ \.ID DISCUSSION 

Jain [8 ' conducted experiments (Tables 3 and 4} using 
brass as material for tool and low alloy steel castings and lov/ 
alloy steel forgings as materials for anode. The shape and 
size of the profiles were measured by preparing a cast replica 
of the hole and projecting it on a tool maker's microscope. The 
method by which anode comer radius was measured, has been des- 
cribed in detail in[^82. Tools have been designed for producing 
profiles obtained experimentally. The designed tool shapes have 
been compared with those experimental tools (Figures 4.5 and 4.8). 

4.6.1. ELECTRnCHB'IICAL DRILLING 

Figure 4.5 shows a comparison between the designed tool 
shape and the tool shape used during experimentation for the 
case of FCD with bare tool. Here, the required and obtained work 
shapes are coinciding because, the design procedure is continued 
till the deviation between the required and obtained work shapes 
lie within an accuracy of 5Hm. 
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Kinks, just opposite to the unevenness in the experi- 
mental work profile, are observed in the designed tool shape. 

This is mostly observed at the end of the transition zone and 
near the top surface of the drilled hole. Kinks obtained during 
designing of the top of the tool (Figure 4.5(b)) is mainly 
because of the curved surface at the top of the work surface 
produced due to stray current attack. 

Figure 4.6 shows the designed bare tool for producing a 
straight sided hole during FCD. It is important to note that 
the designed tool shape may not produce the straight sided hole 
with the specified accuracy, because of the following reason. 

The stray current attack on the work surface is more 
with the bare tools when compared to bit type or coated tools. 

But in the anode shape prediction model, there is no provision 
to account for the stray current attack. 

4.6.?. Bi.FTrnocneacAL bit drilling 

'•igure 4.7 shows the designed tool bit for producing a 
straight sided hole. The design is based on the assumption that 
the voltage remains constant throughout the machining. In other 
words, the designed tool shape will produce a straight hole only 
when, voltage fluctuations during the machining are arrested. 

Figure 4.8 shows a comparison of designed tool bit and 
the one used during experimentation for the case of electroche- 
mical bit drilling. Difference can be observed between the 
obtained and required’ work shapes at the top of the drilled hole 
(Figure 4.8 (a & b)). This is again, because of the same reason 
that the stray current attack is not considered in anode shape 

prediction model. 
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4.6.3. nLECTR0CHB.!l3AL BIT BORING 

Figure 4.9 shows the designed tool bit for enlarging 
and straightening the predillled hole in case of ECBB. As the 
experimental data are not available, the comparison between the 
designed and experimental tools is not made. 

4.6.4. GENERU BTSONGSIOM 

The discrepancy between the designed tool profile and 
the tool profile used during experimentation is attributed to 
the following reasons: - 

(1) The main reason is the inaccuracy of anode shape predic- 
tion model. 

(2) The discrepancy in the transition zone is mainly due to 
the inaccurate measurement of overcut in this zone. The 
Q^ercut at the end of transition zone was measured by 
replica technique. Taking overcut at the end of transition 
zone as radius, an arc was drawn in the whole transition 
zone, assuniing that the overcut in transition zone remains 
constant. But in actual practice, the lEG at the beginning 
of the transition zone is normally smaller than that at 

the end of transition zone* 

•Is discussed earlier, the accuracy of the tool design 
model mainly depends on the accuracy of anode shape prediction 
model. The factors responsible for the poor accuracy of the 
anode shape prediction model (SBFET22) are as follows* 

(1) Stray current attack is not considered. 

(2) Electrolyte flow velocity distribution in ISG is obtained 
from the continuity equation alone. 
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(3) Correct valiies of valency of HC dissolution, machining 

efficiency, accurate electrochemical equivalent of alloys 
and behaviour of gas hubbies distribution in I EG etc* are 
not known* 

4.7. C!';XCGIO:iE 

There is a good agreenent between the designed tool 
shape and the experimental tool shape for the same machining 
conditions. However, the accuracy of the proposed correction 
method for tool design in ECf’, depends mainly on the accuracy 
of anode shape prediction model. The designed tool shape is 
not the complementary shape of the desired work shape in all 
the cases, i.e., ECC, u-gp and ECBB. 
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CHAPTER V 


SCOPE FOR FUTURE UORIC 


The work embodied in this thesis can be extended in 

the following directions; 

(1) An attempt should be made to improve the accuracy of the 
designing model by improving the anode shape prediction 
model in the follovidng ways, 

(a) The effect of the stray current attack on the anode 
profile should be considered. 

(b) flore accurate electrolyte flow velocity distribution 
in IHG should be obtained by solving Navier-Stroke’ s 
equations. 

(2) For the practical applications (for example, elliptical, 
parabolic, shape drilling, etc.) the tool design models 
based on three dimensional analysis should be evolved. 

(3) This tool design procedure can be extended to external 
shaping as well. 

(4) Desired work shape can be obtained fram the given tool 
shape and machining conditions provided independent para- 
meter like voltage or feed rate or both can be accurately 
controlled during the process. For the purpose of real 
time control of these parameters, microcomputer capabi- 
lities can be utilized. The beauty of the proposed method 
( \ppendix-XI) is to use the sane tool for producing diff- 
erent shapes and sizes of workpieces. 
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Such analytical models can also be applied for the 
design of segmented tools (Figure A. 4) whose shape and size 
‘ will be changing with time according to the designed program. 
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APPEMDIX I 
CCT.1PUTATIONAL DETAILS 


A.1 HEAT THANSFER 

The flow chart for the model TFET22 is shown in Figure 
3.5* The other computational details are as follows. 


HEAT GEHERATIOM (fl) 


For a rectangular element of length *dx’, height ’y’ 
and thickness in lEG (Figure A.l(a)) heat generated due to 
Ohmic heating, 


But 


0 0 o 

H = I R S= V /R = KV^ 

K . 4 K .dx.B 

o o 


( A« 1 ) 
(A.2) 


where ^ = cross sectional area perpendicular to the current 

flow. 

Substituting eq. (^.2) in (A.1), the heat generated per 
unit volume for the element, 


H 


H/y.B.dx 



(A.3} 


It is assumed that the heat generated per unit volume 
for all the triangular elements inside the rectangular element 
remains the same. 

As the machining progresses, the work surface becomes 
tapered (Figure A.l(b)). In that case, the element *EFGH' is 
approximated as rectangular element of height y which is given by. 


y 


(y-i + Y2)/2 


(A.4) 
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Again for all triangular elements inside the element 
’ EFGH' , H is assumed to be the same, which is given by eq« (A.3)» 

DEIEF DESCRIPTION OF SUBROUTINES USED 

(I) SUBROUTINE DATIN: 

This subroutine reads in and prints out the input data, 

(II) SUBROUTINE ABC? 

This applies boundary conditions (specified voltages). 

(III) SUBROUTINE BANSOLt 

This solves the sparse, banded, symmetric matrix and 
gives the field potential distribution within the lEG. 

(IV) SUBROUTINE ELBAXTs 

This generates elemental conductivity, mass and heat 
generation matrices. 

(V) SUBROUTINE ASET'.'.BI.J 

This assembles elemental matrices into global form. 

(VI) SUBROTJTINE 3ELf*AT: 

This calculates boundary matrices. 

(VII) SUB RO UT IN E BASMBL : 

This assembles element boundary matrices into a global 

form. 

(VIII) SUBROUTINE rHTf'.ULs 

This performs some matrix operations required for solving 
simultaneous linear, differential equations. 

(IX) SUBROUTINE TBOWD: 

This applies temperature boundary conditions. 

(X) SUBROUTINE F04RF: 

This is a system subroutine used for solving system of 
simultaneous linear equations. 



82 


(XI) SUBROUTINE CDTC: 

This calculates current density (J), conductivity (K) 
at all the nodes, 

(XII) SUBRDTJTINH INET.,GPJ 

This calculates lEG for the case of zero and finite feed 
rate conditions for a node at any instant of time ’t’, 

(XIII) SUB ROUT INE B\VD : 

This determines value of bandwidth. 

(XIV) SUBROUTINE ZERs 

A subroutine that initialises the matrix. 

A, 2. TOOL DESIGN 

Figure 4.4(e) shows the flow chart of the model DES22. 
Some computational details have been presented below. 

’Ihile modifying the tool shape in the side zone, the 
correction has been applied in a direction perpendicular to the 
tool axis (.along y. Figure A, 2). But the correction at any node 
should be applied in a direction normal to the woric surface at 
that node (along Figure A, 2), because the machining takes 

place in that direction. Further, any point ’p* (on the work 
surface), in the beginning of a computational cycle moves to 'p'' 
in a normal direction after that computational cycle (Figure A. 2). 
The movement in the x-direction ’pQ’ (Figure A.2) has been 
neglected in the anode shape prediction model. 

From the Figure A. 2, 

PP' = I f cos9| At 


(A.5) 
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I TC 

X movement, pO *= gr-— - f cos© At cos© (A. 6) 

? V'l ^ 

y movement, p'Q = f cos©] At sin© (A.7) 

L * w 

Similarly, the point *M* on the woik surface moves to 
Fi' (Figure A. 2) and so on. The line joining all such points 
like p' , m’ , .... will give new work surface, after a computa- 
tional cycle. 

If a node ’o' on the tool moves to ’o’’ in a computa- 
tional cycle time of At sec. the new noirmal gap corresponds to 
o'f y^^, is givai by 

y^^ ^ni 1. ■'? P ~ ^ cos©^ A t {A.8} 

where y^^^ = initial normal gap (corresponds to *o*). 

But, in the present model the point 'p' is assxmed to 
move to ’p”' (and f' to *■”* etc.) in the y-direction (Figure A, 2). 
The new lEG at node o’” in the y-direction, y^ (p’’o’” ), is 
calculated as, 

Yo = Yi f cosej At (A.9) 

where = initial lEQ (o'”p) in the y-direction. 

However, for very small computational cycle times 
(30 sec.), the tool movement (oo* ) will be very small. There- 
fore, the x-movement (eq. (A.6)) can be neglected. In this case 
the error incorporated by eq. (A.9) reduces to ’f cos© t' . 

Even though this error is very small (for f = 0.006 mm/s. At - 
30 s, © s= 88^*, error = 0.002 mm), this goes on increasing as it 
is of repetitive nature. 
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For the clearence type of the predrilled hole, the 
front and transition zones have been neglected. The finite 
element discretization starts from the second computational 
cycle onwards as shown in Figure A. 3. In this case, the current 
density for the first computational cycle is calculated using 
the following equation: 

J = K(V - AV)/r^ In(^) 

t 

where r « radius of the work 
w 

r^ = radius of the tool* 

In the anode shape prediction for the full interference type 
(Figure 4.4(b)) of predrilled holes, all the three zones are 
considered in the same way as in case of FCBD. ^or the partial 
interference type of predrilled holes (-igure 4.4(c)), the 
transition zone is considered from the interference point 
onwards. 

For all the cases, the electrolyte flow velocity distri- 
bution is obtained by the continuity equation i.e. U = Q/A. 

But this is not accurate for the case of ECB. 




(a) 



A i AS ^.1 


{ iC’ i\A iAlJ Ut.AI SO?'^ . 
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APPENDIX II 
TOOL DESIGN 


Conventional EC?.! tool design starts with selecting opti- 
mum machining conditions. Next step is to predict tool shape 
for producing the required work shape, while machining under 
prescribed (or optimum) machining conditions. In the subsequent 
stages, the other design aspects like electrolyte flow path, 
insulation etc., are considered. 

The tooling design problem can be attempted from the 
other direction also i.e. controlling the machining conditions 
for producing required work profile, while machining with a 
fixed tool geometry. Apart from the shape of the tool, the work 
profile mostly depends on the independent parameters like vol- 
tage, feed rate and machining time. For a fixed tool geometry, 
the desired work shape can be obtained by varying the voltage 
or feed rate. In case of voltage, the desired work shape can 
be obtained, 

(1) by varying the voltage with time 

(2) by varying the voltage along the tool. 

Voltage can be varied easily with time by suitable means 
viz., by the use of microprocessor or adaptive controls. Figure 
A.5 shows the different analytical work profiles for different 
voltage variations (Figure A.5). The effect of voltage variation 
with time on work profile can be clearly seen from the Figure A.5. 
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Varying the voltage along the tool in a required 
fashion can be possible with the microprocessor and adaptive 
controls. In this case, the tool has to be made of number of 
small bits separated by very thin laminations of insulating 
material (Figure A.4). 
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TABLE j 1 

Thermal Properties of Tool, 'Vork and Electrolyte 


— * 

f 

Material } 

f 

1 

f 

t 

1 

I 

1 

t 

K 

’'J/cm-''C 

■~i 

1 Density (f) 

I 

; gm/cm^ 

t 

f 

t 

- , ^ 

“I 

{Specific heat (C) 

I 

i ^7S/gm^^C 

f 

t 

* ■ r - - 

Copper 

3.8 

8.954 

0.3831 

Steel 

0.6 

7.81 

0.46 

Electrolyte 

0 # 06 

1.0 

4.18 
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TABLE s 2 

Experimental Conditions and Temperature Distribution during Ea/. 
with Plane Parallel Electrodes :;22j 

Total machining time = 300 sec., = 0.1412 

C ss Copper, S s= Stainless steel 


Job 

No. 

1, 

V 

(Volts) 

f 

(mm/ s ) 

Q 

(lit/min) 

Y 

o 

(mm) 

Temperature at distance (<'C) 

X«2.5 =6.5 =17.5 =21.5 ram 

C-1 

12.0 

0.64 

1.5 

0.60 

37.2 

39.3 

42.0 

43.0 

s-i 

12.0 

0.64 

1.5 

0.60 

32.3 

34.7 

40.9 

42.0 

02 

16.0 

0.64 

0.9 

0.82 

36.2 

36.3 

39.8 

41.0 

S-2 

16.0 

0.64 

0.9 

0.82 

32.3 

33.6 

39.5 

40.8 

C«**3 

20.0 

0.64 

2.1 

1.10 

34.2 

36.0 

38.3 

40.8 

S-3 

20.0 

0.64 

2.1 

1.10 

32.3 

34.7 

38.3 

40.8 

C»»4 

12.0 

0.88 

1.4 

0,50 

37.3 

39.8 

42.1 

44.0 

S'»»4 

12.0 

0.88 

1.4 

0.50 

31.3 

33.7 

— 

36.7 


16.0 

0.88 

1.5 

0.60 

37.3 

39.8 

45.6 

48.0 

S-5 

16.0 

0.88 

1 .5 

0.60 

32.5 

35.0 


45.6 

C-6 

20.0 

0.88 

1.7 

0.75 

37.3 

39.8 

44.0 

46.9 

S«“6 

20.0 

0.88 

1.7 

0.75 

32.3 

34.7 

42.0 

45.0 

07 

12.0 

1 .13 

1.3 

0.40 

42.1 

44.0 

49.1 

50.4 

S-7 

12.0 

1.13 

1.3 

0.40 

34.8 

•• 

46.6 

52.5 . 

C-8 

16.0 

1.13 

1 .4 

0.50 

46.6 

50.1 

53.6 

57.1 

3-8 

15.0 

1.13 

1.4 

0.50 

34.8 

39.5 

— 

53.6 

09 

20.0 

1 .22 

1.6 

0.65 

46.7 

47.9 

53.7 

57.1 

S-9 

20.0 

1 .22 

1 .6 

0.65 

36.0 

39.5 

52.0 

56.0 

O10 

12.0 

1 .22 

3.0 

0.60 

30.3 

31.5 

32.8 

35.2 

3-10 

16.0 

1 .22 

1.4 

0.50 

36.2 

39.8 

52.6 

58.4 

■oi 1 

16.0 

0.64 

3.2 

0.84 

31.5 

mm 

35.2 

36.3 

S-11 

12.0 

1 .22 

1.2 

0.40 

33.8 

- 

46,3 

51.5 

012 

20.0 

0.64 

3.1 

1.10 

32.5 

33.8 

36.2 

39.9 

S-12 

20.0 

1.22 

1.6 

0.60 

38.3 

40.8 

53.6 

61,7 
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TABLE J 3 

Experimental Conditions used during ECD of 
Cast Low Alloy Steel rs'j 

f ss 0.00571 rnm/s, s= 1.0 mm 


Work pi ece 
No. 

(ohm~^ nmn'"'* ^ 

E 

V 

(Volts) 

I 

c 

(,'\mperes) 

^ c/ 

Tool-1, r^ 

n 3.79 mm, r^^ = 

1 ,246 mm. 

u = 5.69 m/s 


201 

0.00790 

14.348 

12.52 

83,70 

413 

0.00725 

17.75 

13.1 

90.15 

408 

0.00725 

15.31 

16,24 

91.70 

409 

0.00725 . 

21.84 

14.14 

95.00 

406 

0.00725 ' 

1 2.94 

13.28 

86.9 

Tool— 3, r^ 

c 4.54 mra, r^^ s= 

1 .57 mm. 

u = 5.263 m/s 

- 

407 

0.0063 

12.93 

15,0 

93.16 

421 

0,0070 

14*89 

16.1 

96.20 

424 

0.0060 

18.5 

19.86 

92.0 

425 

0.0066 

20.49 

18.80 

9l .86 


Tool-4, 

r^ = 5.03 mm, r^^ 

= 2.1 3 mm. 

u = 5.37 m/s 


429 

0.0060 

11.57 

18.53 

94.0 

401 

0.00547 

12.975 

22.15 

91.82 

418 

0.00532 

15.95 

18.4 

92.5 

427 

0.00532 

17.62 

- 

* 

5 402 

0.00600 

15.87 

18.4 

93.7 
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TABLE t 4 


^•achining Conditions used during EGBD and SCBB of Low Alloy 
Steel Castings and Low Alloy Steel Forgings *‘8^' 


iS. ss 
0 

0.0007 

-1 -1 

mm ’ , 

^0 = 

1.0 mm, 

Q »s 

5.37 litres/min 

Job 

No. 

V 

(Volts) 

r 

c 

(Amp. ) 

f 

(xlO"^), 

mm/s 

^t 

(ram) 

tc 

(mm) 

h. 

b 

(mm) 

Total 
machining 
time 
(sec. ) 

627 

3.88 

25.67 

3.70 

8.9 

0.5 

5.0 

2400 

635 

7.36 

17.70 

4.50 

6.0 

0 • 5 

4# 0 

2400 

637 

7.64 

17.17 

6.19 

5.5 

1.5 

3.8 

2400 

721 

6.18 

17.02 

5.33 

6.0 

0,5 

3.0 

2400 

725 

7.55 

13.80 

3.70 

6.00 

0.5 

3.0 

2400 

733 

7.60 

16.70 

4.50 

6.0 

0.5 

6.0 

2400 

735 

6.75 

19.41 

5.33 

6.0 

0.5 

5.0 

2400 

736 

6.47 

15.75 

4.50 

6.0 

0.5 

2.0 

2400 

737 

7.33 

16.50 

4.50 

6.0 

2.3 

5.2 

2160 

406 

7.30 

13.66 

4.50 

5.5 

0.5 

2.0 

2400 

413 

8.07 

12.79 

4.50 

5.5 

0.5 

2.0 

1800 

414 

6.30 

15.77 

5.33 

5.5 

0.5 

3.0 

2400 

424 

6.77 

15.07 

6.18 

5.5 

0.5 

3.0 

2400 


* - Inclusive over potential. 
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TABLE * 5 

Experimental Observations About Overcut Inside 
Zone during ECO C8j 


1 — 

Job No* { 

1 



Overcut, 

mm 



Nepth 1 

X, mm 1 

413 

407 

421 

425 

418 

427 

, 

0,0 

2,49 

1.79 

1.71 

2.26 

2.34 

2.32 

0.5 

- 


1.54 

1.76 

- 

1.73 

1.0 

1.9 

1.37 

1 .56 

1.695 

1,93 

1.65 

1.5 

1>r 

- 

1.4 

1 ,61 

- 

1 .63 

2.0 

1.61 

1.26 

1.45 

1 .6 

1.64 

1.44 

2.5 

am 

- 

1.31 

1 .51 

- 

1 .51 

3.0 

1.42 

1.17 

1.36 

1,52 

1.475 

1.34 

3.5 

- 

- 

1.13 


- 

1.42 

4.0 

1.27 

1.05 

1,25 

1 .43 

1 .36 

1.12 

4.5 

- 

- 

1.02 

1,33 

- 

1.31 

5.0 

1.16 

0,94 

- 

1 .34 

1.25 

1.02 

5.5 


- 

0.91 

1 .23 

- 

1.18 

6.0 

1.1 

0.84 

0.92 

1.11 

1.13 

0.91 

6.5 

- 

- 

0.71 

1.1 

- 

1.03 

7.0 

0.96 

0.66 

0,715 

0.91 

0,98 

0.61 

7.5 

8.0 

8.5 

0,35 

0,62 


0.60 

0.84 

0.725 

0,61 
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TABL’H s 6 

Experimental Observations about Overcut Inside Zone 
during EGBn and ECBB',*8'; 


r~ 

Job No . J 

627 

637 

707 

406 

413 

^epth } 

in mm } 




Overcut, ram 



0,0 

1 ,060 

1.335 

0.87 

2.00 

2.01 

0,5 

1.015 

0.825 

0.655 


- 

1,0 

0.895 

- 

0.785 

1 .59 

1.41 

1.5 

0,980 

0.675 

0.655 

- 


2,0 

0.825 


0.78 

1.36 

1.2 

2.5 

0.945 

0.670 

0.645 

, - 

- 

3,0 

0.735 


0.765 

1.21 

0.95 

3.5 

0,780 

0.6725 

0.655 

- 

mm 

4.0 

0.650 

- 

0.765 

1.10 

0.75 

4.5 

0.665 

0.665 

0.65 


- 

5.0 

0.565 

- 

0,765 

1.02 

0.66 

5.5 

0.555 

0,65 

0,65 



&.0 

0,440 

- 

0.725 

0.95 

0.55 

5,5 

0,400 

0,64 

0.625 

- 

- 

7.0 

0.350 

mm 

0,61 

0.88 

0.5 

7.5 


0.59 

0,505 

- 


8.0 


0.54 

0.475 

0.79 


8.5 


0.49 


- 


9.0 




0.62 




